Time-optimal rotation of a spin 1/2: application to the NV center spin in diamond 
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We study the applicability of the time optimal bang-bang control designed for spin-1/2 [U. Boscain 
and P. Mason, J. Math. Phys. 47, 062101 (2006)] to the rotation of the electron spin of a nitrogen- 
vacancy (NV) center in diamond. The spin of the NV center is a three-level system, with two levels 
forming a relevant qubit subspace where the time- varying magnetic control field performs rotation, 
and the third level being idle. We find that the bang-bang control protocol decreases the rotation 
time by 20-25% in comparison with the traditional oscillating sinusoidal driving. We also find that 
for most values of the bias field the leakage to the idle level is very small, so that the NV center 
is a suitable testbed for experimental study of the time-optimal protocols. For some special values 
of the bias field, however, the unwanted leakage is greatly increased. We demonstrate that this is 
caused by the resonance with higher-order Fourier harmonics of the bang-bang driving field. 

PACS numbers: 02.30.Yy, 76. 30. Mi, 03.67.-a 



I. INTRODUCTION 

Studying the dynamics of quantum spins helps to bet- 
ter understand many fundamental issues of quantum 
physics, and is important for applications in quantum 
information processing,— coherent spintronics,^ and high- 
precision metrologyj^ Recently, spins of the nitrogen- 
vacancy (NV) centers in diamond have emerged as 
a particularly strong candidate for long-range quan- 
tum communication,^ scalable quantum computation,'^'^ 
and high-sensitivity magnetometry with nanoscale 
resolution. NV centers also present an excellent plat- 
form for exploring such fundamental problems as dynam- 
ics of quantum spins coupled to their environment f^^rJ^ 
quantum control and dynamical decoupling of single 
solid-state spinsJ^"— This system has become very pop- 
ular due to a favorable combination of properties: a 
NV electronic spin has long coherence time^^'iiiii'i^ a 
single NV center can be optically initialized and read 
out,^^'^° and can be manipulated both optically^-'^ and 
magnetically. •^S 

In this paper we consider fast manipulation of the spin 
of a single NV center by externally applied magnetic field, 
which is important for many prospective applications. 
The magnetic field can be quickly varied in time with 
modern electronics, but its magnitude is always limited, 
so the problem of time optimal spin control naturally 
arises: how to shape the time profile of the control field so 
that the desired spin evolution is accomplished within the 
shortest time. Time-optimal control is a well developed 
science nowadaySf^l and the optimal control of quantum 
spins has been extensively studied in various contexts, 
from general mathematical considerations of the opti- 



mal control of two-level systems 



24-29 



to numerous ap- 



plications, e.g. in the areas of magnetic resonance) ^^i^° 
and quantum information processing4ir"— In most cases, 
when the system under consideration is sufficiently com- 
plex, one has to use numerical methods to find the op- 
timal controli ^"'^^ There are, however, rare cases where 
the exact solution is known, and its optimality is rig- 



orously proven. ^^'^^ It is of interest to see whether NV 
centers in diamond may present a suitable system for ex- 
perimental test of such fundamental results of the control 
theory. Besides basic interests, such protocols can speed 
up control of NV spins in future applications. 

We consider the exact solution for time-optimal control 
of a spin 1/2, placed in a static magnetic field along the 
2;-axis, with the corresponding energy splitting 2E be- 
tween the states | t) and 14-)- As in standard resonance 
experiments,'^^ a driving field Bx{t) of limited amplitude, 
< M, can be applied along the a;-axis to rotate 
the spin from | J,) to | f)- The traditional solution is 
to use driving which oscillates with the spin's Larmor 
frequency 2E/h: for small driving, M <^ E, within the 
rotating wave approximation (RWA), it implements the 
desired rotation. However, this approach is not optimal; 
also, at strong driving the RWA breaks down, making 
the spin dynamics intricatej ^^''^^!'^^ It has been provenSl 
that for a single-axis driving the bang-bang control is op- 
timal, with the driving field Bx{t) switching between its 
extremal values. ±M 

Here we investigate whether this exact solution is ap- 
plicable to the spin of a NV center. A NV center has 
spin S = \ with the anisotropy splitting D = 2.88 GHz 
between the levels ms — ±1 and mg = 0. Static exter- 
nal field applied along the anisotropy axis further splits 
the levels mg = 4-1 and ms = — 1. In most applica- 
tions, the states mg — and ms = — 1 form the relevant 
two-level system (qubit, or a pseudo-spin 1/2) while the 
third level ms = 1 remains idle. However, when the op- 
timal control designed for a spin-1/2 is applied to the 
NV spin S = 1, the idle third level can become occupied. 
Protocols minimizing the leakage to higher levels have 
been developed before, using numerical tools, for differ- 
ent physical systems and situations 42r— But numerical 
optimization does not guarantee global optimality, and is 
model-dependent. Thus, it is interesting to consider the 
opposite approach, and study applicability of the known, 
rigorously optimal spin-1/2 protocol, to a NV spin rota- 
tion (although such protocols could be sub-optimal when 
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FIG. 1: (a) Energy levels of the electronic spin S = 1 of a NV 
center. The states ms ~ and mg = — 1 form the relevant 
qubit subspace. (b) The position of the energy levels when 
the high order harmonics in the control field can induce the 
resonant leakage to the idle level \ms = 1). 



the full 3-level system is considered). We show that for 
most regions of the relevant parameter space, the leak- 
age to the idle level is negligible, so that NV center does 
present a good testbed for the spin-1/2 protocol. This 
also means that the known time optimal protocol can be 
safely used to manipulate the NV center's (pseudo)spin, 
and can noticeably, by ~ 25%, speed up its rotation com- 
pared to RWA. Moreover, we investigate the cases where 
the known solution for spin 1/2 is not applicable, and 
explain the reasons for the increased leakage to the idle 
NV spin state. 

The rest of the paper is organized as follows. In Sec. II, 
we describe in more detail the spin Hamiltonian of the 
NV center, the time optimal spin-1/2 protocol, and our 
results. Discussion and conclusions are given in Sec. III. 



II. RESULTS: BANG-BANG CONTROL OF THE 
NV CENTER AND ITS PERFORMANCE 

The nitrogen- vacancy color center in diamond consists 
of a substitutional nitrogen atom with an adjacent va- 
cancy. Its many-electron orbital ground state has a to- 
tal spin^"^ S = I. The combined effect of the spin- 
orbit and spin-spin interactions leads to a single-axis spin 
anisotropy, with the anisotropy axis directed along the 
(111) crystallographic direction (symmetry axis of the 
NV center) which we take as the ^-axis. As a result, the 
states ms = ±1 and ms = are split hy D = 2.88 GHz 
(everywhere in this paper, we assume h = 1). The exter- 
nal static field applied along the z-axis further splits the 
states ms = 1 and ms — —1, as shown schematically in 
Figure [T] (a), and the total Hamiltonian of the system is 



Ha = DSl + BoS, + B,{t)S., 



(1) 



where Bx{t) is the control field to be optimized. When 
NV spin is used, for instance, as a qubit for quantum 
information processing, only the two levels ms = and 



ms = — 1 are employed. The control field B.j;{t), whose 
amplitude is bounded \Bj;{t)\ < AI, is applied to drive 
the spin from the initial state ms = to the final state 

ms = -1. 

Within the relevant subspace spanned by these two 
eigenstates, the system is equivalent to a (pseudo)spin 
s = 1/2. The Hamiltonian for this spin, omitting the 
irrelevant energy shift, is 



H ^ Es'' + V2B^{t)s'' 



(2) 



where E — {D + Bo)/2, and Bx{t) should steer the spin 
from I I) to I t)- The optimal control for this prob- 
lem has been thoroughly investigated in Ref. [27]. The 
time-optimal rotation of the (pseudo)spin is implemented 
by the bang-bang control, where the control field B.j;{t) 
switches between +M and —M. The term bang here, 
following the standard terminology, denotes the interval 
between two neighboring switchings, i.e. during a bang 
the control field is constant. When the maximal control 
field is strong, M > V2E, the time-optimal rotation is 
achieved with exactly one switching^? 



M 



0<t<ti 



'M ti<t<T 



where the switching time ti — 
SiTCCos{2E^/AP)]/{2y/E^ + and the 

control time is 



T = 



^E^ + AP/2 



(3) 

[tt ± 
optimal 

(4) 



Corresponding spin trajectory on the Bloch sphere is 
rather simple: the spin vector starts from the south 
pole of the sphere (| \) state), rotates about the axis 
2Ez + \/2Mx until it hits the equator, then continues to 
rotate about the axis 2Ez — \f2Mx all the way to the 
north pole (| f) state). Here x and z are the unit vectors 
along X and z axes, respectively. Note that the point 
i? = is singular: a static field of maximum amplitude 
would perform the necessary rotation twice faster than 
the protocol of Ref. [13. However, this can be used only in 
very special situation, in the nearest vicinity (E <C M) 
of the point _E = 0, so we do not consider this singular 
case below. 

For weak driving M < ^/2E, the duration of all bangs 
is the same, except for the initial and the final bangs: 

B^it) = {-l)^-^M, te[U + {k - l)tra,U + ktm) (5) 

where i^, tm, and tf denote the durations of the initial, 
middle, and the final bangs, respectively, nswt is the 
number of times the field switchings between -\-M and 
— M, and fc = 1, 2, • • • riswt + 1 denotes the bangs. There 
are also theorems relating tm and tf toU, and limitations 
on the number of the field switchingsi^ 

Instead of the protocol of Ref. [23, a traditional ap- 
proach for spin rotation employs a driving field along 
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FIG. 2: Spin rotation time as a function of the static 
field Bo for different amplitudes of the control field, M — 
50, 200, 400, 600 MHz from top to bottom, respectively. The 
rotation time at the peak of each curve is (accidentally) equal 
to tRWA, and is denoted by the points of different shapes. The 
number on the right of each curve indicates the ratio of the 
control time to tawA at large bias Bo- 



the X axis oscillating at Larmor frequency, Bx{t) = 
M cos{2Et). This control works well in the weak-driving 
regime, M <^ E. Within the rotating-wave approxima- 
tion (RWA))^ when only the secular terms are retained 
in the control field, the oscillating field is equivalent to a 
driving field with magnitude M/2 which rotates in the x- 
y plane in step with the Larmor precession of the spin s. 
The corresponding rotation time is ^rwa = v^tt/M. In a 
general situation, when the control field can be directed 
anywhere in the x-y plane, the rotating field is known 
to provide time-optimal control. However, in reality it 
is very diflacult to create two perpendicular control fields 
with precisely locked directions, amplitudes, and phases. 
Therefore, in the case of the oscillating driving along the 
X axis, half of the control field approximates the opti- 
mal two-axis rotating-wave solution, while the other half 
(corresponding to the omitted non-secular terms in the 
RWA) is wasted. The solution of Ref. '2? utilizes the full 
magnitude of the single- axis control in an optimal way. 
At the same time, for M <^ E, its approximately peri- 
odic structure resembles the oscillating driving, but the 
period (duration of the bangs) is slightly different, and 
its time shape is rectangular rather than sinusoidal. 

In a typical experiment, the maximum magnitude of 
the driving field M is fixed, limited by the setup de- 
tails, while the static external field Bq is adjusted to 
provide the convenient conditions for driving the rota- 
tion of the NV center spin from ms = to ms — —1. 
Fig. m shows the rotation time for the protocol of Ref. [131 
as a function of the bias Bq for different driving ampli- 
tudes M = 200, 400, 600 MHz. Overah, as expected, 
the rotation time decreases approximately as I/A/. At 
Bo = 1028 G the distance 2E = D + Bq between ms = 



and 7715 = — 1 is zero (level crossing point), and the ro- 
tation time is maximal. The corresponding control time 
2tt/ {^/2M) is, accidentally, equal to t^wA, although RWA 
is not applicable here. Farther from the level crossing 
point, the rotation time decreases. At very large dis- 
tance 2E the rotation time almost saturates at the value 
not very far frontal t'ssNA/^/2. Therefore, the protocol of 
Ref. [23 can decrease the rotation time by about 25% in 
comparison with the standard oscillatory Bx{t). 
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FIG. 3: (Color online). Dependence of the final population 
P_i (red line) and Pi (blue line) on the static field Bo, after 
the application of the optimal control for the corresponding 
two-level system, (a) The results for M = 50 MIfz. The inset 
shows the values ln[l — P_i] (red) and ln[Pi] (blue) vs. Bq. (b) 
The results for M = 200 MHz. The inset shows the Fourier 
transform of the bang-bang control, Bx{v), for static bias 
Bo — D/2. The first peak at 1.435 GHz corresponds to the 
desired transition between levels ms = and ms = — 1, while 
the second peak at 4.305 GHz drives the unwanted transition 
to ms ~ I state. 

Fig. El shows the fidelity of the spin-1/2 protocol when 
applied to a NV spin S — 1. Except for a few sharp 
dips at certain bias fields, the final population P_i of the 
targeted level ms = — I stays high, and far from the level 
crossing saturates at the values 0.999, 0.994, and 0.986 
for M = 50, 200^nd 600 MHz, respectively. Thus, the 
protocol of Ref. [23 provides good fidelity for S" = 1 spin 
of NV center for most values of Bq . The loss of fidelity, 
including the neighborhood of the dips, is predominantly 
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due to leakage into the idle state mg = 1, as seen in 
Fig. [21 and the dips in P_i are accompanied by the peaks 
of Fi. 

Such leakage is particularly prominent when the energy 
levels satisfy the condition 



e+1,0 = (2n + l)|e-i,o| 



(6) 



where e±i,o are the energy splitting between ms — ±1 
and ms — 0, respectively, and n = 0, 1, 2 ... In terms of 
Bq, these conditions are 



Bo 



n + 1 



-D or 



1 



-D 



(7) 



The increased leakage is determined by the (approxi- 
mately) square- wave time profile of Bx{t), with (approx- 
imate) period 2<m, see Eq. [S] The Fourier expansion of 
the square wave contains all frequencies (2n + 1)tt /tm 
with integer n, and these features are prominent in the 
Fourier transform of the bang-bang control Bx{t) (inset 
of Fig. 13]) . The harmonics of the driving field come in res- 
onance with the transition between ms = and mg = 1, 
see Fig. [H thus drastically increasing the leakage. To 
demonstrate this, we set the external field at the value 
Bq — D/2 {n ~ 1) corresponding to the middle of the 
pronounced dip in Fig. [31 Instead of the actual bang-bang 
control Bx{t), which leads to P_i = 0.874, we used its fil- 
tered version Bf(t), by cutting off the Fourier harmonics 
with frequencies above certain value Vc- Also, to make 
the comparison fair, the control time for filtered control 
was kept the same as the original bang-bang control time 
(by turning off the filtered control Bf{t) at longer times, 
see Fig. [3]). Fig. [H shows the dependence of P_i on the 
cutoff value: as soon as Vc becomes smaller than 3D/2, 
the value of P_i demonstrates a clear jump, since the un- 
wanted resonant harmonic in the driving field is filtered 
out almost completely, and the leakage to ms = 1 state 
is suppressed. 



III. CONCLUSIONS AND DISCUSSIONS 

NV centers in diamond represent a very suitable sys- 
tem for fundamental studies of quantum spin dynamics, 
and for many quantum spin-based applications. In many 
applications so far, the spin of a NV center was used to 
represent a qubit, with the states ms = and ms = — 1 
forming the relevant two-level subspace, and the state 
ms = 1 remaining idle. We investigated whether the spin 
S" = 1 of a NV center constitutes also a suitable testbed 
for studying the time-optimal control protocol designed 
for two- level systems. We found that by applying the op- 
timal control protocol developed in Ref. 21r, the rotation 
time from ms = to ms — —1 state can be reduced by 
^ 25% in comparison with the standard sinusoidal oscil- 
latory driving (within RWA approximation). We checked 
how important is the unwanted leakage to the idle state 
ms = 1. We found that this leakage, for most values 
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FIG. 4: (Color online). Final population P_i of the ms = — 1 
state as a function of the cutoff frequency Vc, for M = 200 
MHz and Bo ^ D/2 = 1.435 GHz (where the energy gap 
between ms ~ and ms = —1 equals Bo). When Vc becomes 
smaller than e+i,o, the leakage to ms = +1 is suppressed, 
and P_i abruptly increases. The inset shows the time profile 
of the original bang-bang (black line) and filtered (blue line) 
controls. The bang-bang control has ti — tf = 0.1266 ns, 
tm = 0.3497 ns, with nswt ~ 8 switchings. 



of the external field, is negligible, so that the NV center 
indeed could be a suitable platform for studying the time- 
optimal control problems. We found, however, that near 
certain values of the bias field Bq, the leakage to the idle 
state is drastically increased, and we demonstrated that 
this is induced by higher-order harmonics in the bang- 
bang driving. We show that by frequency filtering, this 
leakage can be noticeably suppressed. 

Experimental implementation of the protocol studied 
above on NV centers seems reasonable. Currently, the 
magnetic fields which drive the desired rotation can be 
controlled on timescales of tens of picoseconds (or even 
faster with specialized techniques) . The driving field am- 
plitudes of hundreds of MHz have also been implemented 
in many experiments J^^— Moreover, for faster rotation, 
it is important to mover far away from the level crossing 
point at Bq = 1028 G, so that negative external bias or 
very large external field is more favorable for the rotation 
speedup. 

It is important to note that the focus of this work is to 
check whether a NV spin is a suitable testbed for exper- 
imental studies of the time-optimal protocols. Although 
the protocol studied above does speed up the NV spin 
rotation, it most probably presents a sub-optimal solu- 
tion for the full system S = 1. However, the detailed 
studies of the time-optimal control of three-level systems 
present a difficult mathematical problem, and are still 
lacking. In this case, using an optimal solution within 
a restricted subspace is a reasonable option, which may 
complement the numerical solution of the time-optimal 
problem. 
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